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ABSTRACT 

r^ I We examine the relationship between the mass and x-ray gas temperature of galaxy clusters 

using data drawn from the literature. Simple theoretical arguments suggest that the mass of a 

3/2 

cluster is related to the x-ray temperature as M oc Tx • Virial theorem mass estimates based 
on cluster galaxy velocity dispersions seem to be accurately described by this scaling with a 
normalization consistent with that predicted by the simulations of Evrard, Metzler, & Navarro 
(1996). X-ray mass estimates which employ spatially resolved temperature profiles also follow 
\r-s I a Tx scaling although with a normalization about 40% lower than that of the fit to the virial 

masses. However, the isothermal /3-model and x-ray surface brightness deprojection masses 
0| , follow a steeper ex X^-8-20 gcaling. The steepness of the isothermal estimates is due to their 

Qv ' implicitly assumed dark matter density profile of p(r) oc r^^ at large radii while observations 

Q> ■ and simulations suggest that clusters follow steeper profiles (e.g., p{r) ex r~^"'). 

Q ■ 1. Introduction 

^ ' The relationship between the mass and x-ray temperature of galaxy clusters is a necessary bridge 

between theoretical Press-Schechter models, which give the mass function (MP) of clusters, and the 



observed x-ray temperature function (TP) (e.g.. Edge et al. 1990; Henry & Arnaud 1991; Henry 1997 



^^ ' Markevitch 1998). Theoretical arguments suggest that the virial mass of a galaxy cluster is simply related 



f ' i 3/2 

H , to its x-ray temperature as M oc Tx . These arguments are supported by simulations which show a tight 



correlation between mass and temperature (e.g., Evrard, Metzler, & Navarro 1996; Bryan & Norman 1998 



Eke, Navarro, fc Frenk 1998 ). This suggests that the M-T relationship may also be a relatively accurate 
and easy way to estimate cluster mass. However, the M-T relationship first needs to be calibrated using 
masses estimated by some other means. 

The oldest method of measuring cluster mass is the virial mass estimate based on dynamical analysis 
of the observed velocity dispersion of the cluster galaxies. The existence of the x-ray emitting ICM of 
galaxy clusters allows an independent mass estimate but requires knowledge or assumptions about both the 
x-ray temperature and surface brightness profiles. More recently, strong and weak gravitational lensing by 
clusters has provided a third independent mass estimate. If clusters are dynamically relaxed and relatively 
unaffected by non-gravitational process, these three methods should give the same results. In this paper, 
we concentrate on virial and x-ray mass estimates. 
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In Section g, we discuss the theoretical basis of the M-T relation and the results of cluster simulations. 
We then compare this relation with those using masses based on galaxy velocity dispersions (Section y), 
x-ray mass estimates of clusters with spatially resolved x-ray temperature profiles (Section 4.1), and 
isothermal x-ray mass estimates (/?-niodel estimates in Section |4.2| and surface brightness deprojection in 
Section 4.3). In Section o, we discuss the results and present conclusions. 



2. Theory and Simulations 

For gas that shock heats on collapse to the virial temperature of the gravitational potential, the average 
x-ray temperature 

Tx oc oc M^j; (1) 



where rvir is the virial radius, the boundary separating the material which is close to hydrostatic equilibrium 
from the matter which is still infalling. The coefhcient of this relationship is a complicated function of 



cosmological model and density profile of the cluster (see e.g., Lilje 1992). However, because the infall 
occurs on a gravitational timescale igrav ex p~^'^ , the virial radius should occur at a fixed value of the 
density contrast, defined as: 

^ W) M{< r) 



Pc{z) ^TTpc{z)r^ 



(2) 



where Pc{z) is the critical density. For an Hq = 1 universe, Jvir = IStt^ w 178 but drops to lower density 



contrasts for lower values of fio (e.g., [Lacey &: Cole 1993| ). Since we do not have a priori knowledge of the 
actual value of (5vir, we scale all of our results to 5 ^ 200, which should contain only virialized material, and 
has been used previously by other authors (e.g., [Carlberg et al. 1996 ; Navarro, Frenk, fc White 199E ). 



Evrard et al. (1996) (hereafter EMN) present a M-T relation which seems to well describe simulated 
clusters in six different cosmological models (two VLq = 0.2 and four JIq = 1 models, see their Table 1 for 
details). EMN assume that M oc Tx and then fit the coefficient of the relationship at various density 
contrasts. For 5 sa 200, the coefficient depends only weakly on Og, with a difference of « 5% between the 
r^o — 0-2 and f^o = 1 models while the difference rises to k, 40% at (5 = 2500. 

From EMN's Equation 9 and fitting the normalization using their Table 5 (excluding the 5 — 2500 
values but including both ilo = 0.2 and f^o = 1 points), the expression for mass as a function of temperature 
and density contrast is: 



M{5, Tx) = (1.81 ± 0.23) X W5 



14r-0.266±0.022 



\ke\f ) 
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Mr. 



(3) 



where h is the Hubble constant in units of 100 km s ^ Mpc ^. Note that Equation || shows that clusters in 
the EMN simulations have dark matter density profiles p{r) oc r~^-^ in accordance with the effective slope 



of the universal density profile of Navarro, Frenk, fc White (1996) in the relevant range of radii. 

Other simulations using different cosmological models and codes generally give normalizations similar 
to Equation ^ to within < 20%. Eke et al. (1998) give masses and gas temperatures at a density contrast 
of 100 for their simulations of an f^o — 0.3 J^a — 0.7 model which are well described by Equation |^. Bryan 
fc Norman 1998 give the M-T at (5 = 250 for simulations using variety of cosmological models. Like EMN 
they find the normalization is fairly insensitive to the model used although their normalizations are about 



20% higher. Lastly, Metzler & Evrard (1997) examine the differences between simulations of galaxy clusters 



with and without galactic winds. The M-T relation in the wind models is slightly steeper due to the energy 
injection, with a power law index ^^ 1.6 and has temperatures about 20% higher at the low mass end. 

A redshift dependence is introduced into the normalization of the M-T relation by the definition of 



density contrast since the critical density is a function of fl and the redshift of formation (see e.g., Lilje 



1992; Eke, Cole, & Frcnk 1996; Voit & Donahue 199S). This should not substantially affect our results as 



the samples considered consist mainly of low redshift (z < 0.1) objects and/or have scatter in the mass 
estimates considerably greater than the effect introduced by the redshift dependence. 

Since we have restricted ourselves to low density contrasts and redshifts, we will not discuss gravitational 
lensing mass estimates. Lensing estimates are usually limited to high density contrasts 6 > 3000 (even for 



weak lensing) and to moderate-to-high redshifts. However, Hjorth, Oukbir, fc Van Kampen (1998) have 



reported good agreement between the EMN relation and their sample of eight lensing clusters. They assume 
M oc Tx^^ and an 17c 
EMN normalization. 



M (xTJ and an 17o = 1 cosmology. Their best fit normalization is 0.82 ± 0.38 (rms dispersion) times the 



3. Virial Theorem Mass Estimates 

Assuming that the galaxies are distributed similarly to the total mass, the virial mass of a cluster 
depends on the line of sight projected velocity dispersion of the galaxies, Up, and the virial radius, r^r'- 
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If the entire system is not included in the observational sample, as is common for galaxy clusters. 
Equation Q overestimates the mass. The usual form of the virial theorem (2U -I- T = 0) should be modified 
to include the surface term (2U + T = 3PV) since the surface pressure reduces the amount of mass needed 
to keep the system in equilibrium (sec pirardi et al. 1998; Carlberg et al. 1997). 



Girardi et al. (1998) (hereafter G98) have derived virial masses for 170 nearby clusters {z < 0.15) 



using data compiled from the literature and the ENACS data set pK^atgert et al. (1998)| . They define the 
virial radius to be r^r — 0.002ap h~^ Mpc where ap is in km s~^ and consider only galaxies within this 
radius in the mass estimation. Their quoted masses are generally smaller than previous estimates by < 40% 
which they attribute to stronger rejection of interlopers and a correction factor of '^ 19% accounting for the 
surface term. 



We have cross-correlated the G98 catalog with two catalogs of ASCA temperatures (Fukazawa 1997; 



Markevitch 1998) to obtain a subsample of 30 clusters with at least 30 redshifts within Vmr- Fukazawa 
and Markevitch excluded the center of the clusters to minimize the effect of cooling flows on the derived 
temperature. In cases of multiple measurements, the temperatures generally agree within their quoted 
errors (usually <^ 10%), and we have preferentially used the Fukazawa temperatures. This sample has few 
low temperature clusters, so we have analyzed the archival ASCA data of three additional G98 clusters (see 
Table 0). Cooling flow effects are not important for these three clusters. Table lists clusters in the final 
sample along with their adopted x-ray temperature (column 2), G98 virial radius (column 3), and virial 
mass (column 4). 

The assumption that r^ir = 0.002ap is quite approximate, and the actual relation between r^r and ap 
depends on the cosmological model. The density contrast of the G98 virial masses (= Myir / {4:/3^^rl^^pc)) 
is generally less than 200 with a mean (and standard deviation) of 97 ± 23. Assuming the dark matter 



density p{r) oc r^^-^ in the outer parts of the clusters, which both the EMN simulations and Girardi data 
(at least the galaxy distribution) seem to follow, we have rescaled their masses to 6 = 200. Effectively this 
is just a change of normalization such that the rescaled masses are smaller than My^^ by about 15% (with a 
standard deviation of about 5%). 

Figure shows the distribution of rescaled virial masses versus temperature for this subsample. A 



power law fit using the BCES bisector method of Akritas fc Bcrshady (1996) , which takes into account 
the errors in both variables and the possibility of intrinsic scatter, gives M oc j^'i'i±o.20 ^^y[ quoted errors 
are la unless otherwise stated), marginally inconsistent with the EMN relation. The seven most severe 
of the outliers in this plot are A119, A754, A1656 (Coma), A2256, A2319, A3558, and A2029. Except 
for A2G29, all are known to contain complex velocity or temperature structure. We have marked these 
six clusters in gray in Figure |^. Removing these six clusters from the fit gives M oc y^-53±0-i3 and a 
normalization consistent with the EMN relation (see Table M. Further outlier removal or permutations 
of the temperatures (i.e. using Markevitch instead of Fukazawa temperatures) does not cause significant 
differences in the fit (e.g., power law index is changed by ±0.05). Given the relatively large scatter, more 
clusters with well measured temperatures, especially cooler/less massive clusters, are needed to further 
constrain the relationship between virial mass and x-ray temperature. 



3.1. The Velocity Dispersion Temperature Relation 

To first order, Myir oc a^ since Myir oc a'^r^ir and ryir oc a. Therefore, the virial mass-temperature 
relation is related to the more extensively discussed velocity dispersion-temperature (tr-Tx) relation (e.g., 



Lubin fc BahcaU 1993| ; [Bird, Mushotzky, fc Metzler 1995| ; |Wu, Fang, fc Xu 1998| and references therein). 
Figure shows the tr-Tx relation for our G98 subsample. Excluding the same six clusters as above as fit 
gives a oc yJ-54±o.03^ This is consistent with the observed Myir oc T, 



1.5±0.1 



If galaxies and gas are both in equilibrium with the cluster potential and gravity is the only source of 
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energy, cr oc Tx (e.g.. Bird et al. 1995). Our results are marginally steeper than this and consistent with 
many previous estimates (e.g., 0.61 ± 0.13 from Bird ct al. 1995). However, our fit is shallower than the fit 
G98 give in their paper (0.62 ± 0.04) possibly due to improved x-ray temperatures since they use Einstein 



and Ginga temperatures from David et al. (1993) and White, Jones, fc Forman (1997) . However, our fit is 
definitely shallower than the 0.67 ± 0.09 found by Wu et al. (1998). While they draw a much larger sample 
(94 clusters) from the literature, their sample is heterogeneous and the quality of their data is unclear. In 
contrast, we are using only G98 derived velocity dispersion with greater than 30 member redshifts and 
precise ASCA temperatures. 



3.2. Scatter in Virial Mass Estimator 

The scatter in the virial mass estimator is expected to be quite large because of shot noise due to 
the finite number of galaxies in a cluster and projection effects due to contamination by background and 
foreground galaxies. The scatter in the observed virial mass-temperature relation is then a combination 
of the dispersion in the virial mass estimator (with respect to the true cluster mass) and any intrinsic 
dispersion in the M-T relation. 

Figure g shows a histogram of the ratio of the rescaled virial masses (M2oo(G98)) to the mass expected 



from the best fit relation (Mfa). The mean (or median) is approximately 1.0 with a standard deviation of 
0.31. Figure || excludes the six outliers that were not fitted, including these clusters increases the standard 
deviation to 0.59. The expected scatter in the virial mass estimator has not been widely reported in the 



literature, but Fernley fc Bhavsar (1984) find that in their simulations of galaxy clusters the ratio of the 
virial mass to true mass is 0.97 ± 0.36 {la standard deviation) after removing contaminating background 
and foreground galaxies. This predicted scatter is close to the observed scatter around the fit and suggests 
the dispersion in the virial mass - temperature relationship is primarily due to the scatter in the virial mass 
estimator. 

This is further supported by the distribution of Al2oo{G98)/Mfit as a function of the number of 
redshifts {uz) used to calculate the virial mass (see Figured which also includes clusters with less than 30 
redshifts). The scatter is about a factor of 2 lower for clusters with n^ > 80. This is not surprising as a 
larger number of redshifts increases the accuracy of the virial mass estimator by decreasing the shot noise. 
Together with the results in Figure 3, this further indicates that the M-T relation must have very small 
intrinsic scatter. 



4. X-ray Mass Estimates 

X-ray mass estimates are based on the assumptions that the ICM is in hydrostatic equilibrium and 
supported solely by thermal pressure. With the further assumption of spherical symmetry, the gas density, 
temperature, and pressure are related to the mass by: 



dr 



Pg 



Pn 



G M{< r) 



jimp 



(5) 



(6) 



where k is Boltzmann's constant, and ^irrip is the mean molecular weight of the gas. The enclosed mass at 
a radius, r is then: 



M{< r) 



G^ra-p 



d log pg{r) d log T^{r) 



d log r d log r 

which depends on both the gas density and temperature profiles. The gas density is usually assumed to 
follow a (3 profile: 

2 



(7) 



PsW =Pff(0) 



l+(- 
rc 
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(8) 



Historically, x-ray detectors have either had good spatial or spectral resolution but not both. Measuring 
the actual temperature profiles of clusters has really only become practical with the advent of ASCA and its 
ability to obtain spatially resolved spectra. The combination of using ROSAT to obtain Pg{r) and ASCA to 
obtain T^{r) can yield cluster masses with unprecedented accuracy. However, estimating the temperature 
profile is complicated by the spatial and energy of the ASCA PSF, and mass estimates have only been 
reported for a few of the best observed clusters. The difficulty of measuring the temperature profiles means 
that there are far larger samples of clusters for which only the average x-ray temperature and isothermal 
mass estimates are available. 



4.1. Mass Estimates with Spatially Resolved Temperature Profiles 

No large catalog of clusters with masses measured using spatially resolved temperature profiles has 
been published. Therefore, we have searched the literature (including conference proceedings) to obtain a 
sample of 12 clusters with masses measured using known temperature profiles. These clusters are presented 
in Table H with the adopted average x-ray temperature and 90% confidence limits (column 2 or table notes), 
largest radius in which the mass was given (column 5), the mass within that radius (column 6), and the 
reference from which the data was taken (column 10). If given, we have used global temperature values 
given by the respective authors. Otherwise, we have taken them from catalogs of ASCA temperatures as in 
Section 0. The formal errors on these mass estimates are small as they are well constrained by the density 
and temperature profiles. However, the systematic uncertainties (i.e. uncertainties in the ASCA PSF and 
effective area) are much more difficult to quantify. For fitting purposes, we have chosen not to weight the 
fit with any mass errors, only with the errors in temperature. 



Several of the clusters used warrant some comments. The masses for A496 and A2199 from Mushotzky 



et al. (1994) were derived without ASCA PSF corrections. However, the corrections are not large for these 



clusters as they are relatively cool and only the central field-of-view of the telescope was used. A2256 is 



known merger, but [Markevitch fc Vikhlinin (1997)| argue that the subcluster is physically well separated 
along the line of sight and has not disturbed the bulk of the primary cluster's gas. Although this cluster 
was considered an outlier in the virial mass fit, this can be attributed due to contamination of the optical 
velocity dispersion measurement. 

As with the virial masses, we have rescaled the masses to a density contrast of 200 assuming the dark 
matter density profile p(r) ex r"^"'. On average, the effect of rescaling is to increase the masses by an 



average of 20% with a standard deviation of w25%. There is some support for using this profile. Markevitch 



et al. (1996)1 and [Markevitch fc Vikhlinin (1997)| report similar profiles for A2163 and A2256. While other 



authors in Table g do not give density profiles, some (e.g., phashi 1997 ; Sarazin, Wise, fc Markevitch 1998 ) 



contain plots of the mass as a function of radius which are consistent with the p(r) ex r~^'^ profile. Lastly, 
the clusters themselves seem to obey this profile. Figure || show the density contrast (effectively the average 
density) as a function of radius. We normalize the radius to EMN's r2oo, the radius at which the density 
contrast is 200, so that we are comparing similar scales in different clusters but this makes little difference. 
An unweighted fit indicates that p{r) ex r^^-"^^^-^. 

In Figure ^, we plot the rescaled masses versus temperature. Excluding A3526 (Centaurus), the best 
fit is M oc x^-48±o.i2 ^^|. ^[^Y^ g^ normalization about 40% lower than that of the EMN relation. Given the 
heterogeneous nature of the sample, the dispersion around this fit is surprisingly small {<^ 10% in mass) 
indicating that the intrinsic correlation between temperature and mass is quite tight. Using a different 
density profile to extrapolate the mass to a density contrast of 200 has the tendency to increase the 
dispersion in this fit (i.e. about 25% for p{r) oc r~^) but does not have much effect on the power law index 
or normalization of the fit. 

The lower normalization than that found using the EMN relation or virial masses may reflect 
systematics in the masses derived using the temperature profiles, or it could be a problem with the 



simulations and systematics in the virial mass determinations. Aceves fc Perea (1998) have found that the 
virial mass estimate can either overestimate or underestimate the mass depending on the aperture of the 
region sampled. However, a direct comparison of the virial mass and x-ray masses for the 9 clusters with 
temperature profiles masses (M2oo(T(r))) and G98 virial masses (M2oo(G98)) (Figure 0) shows no clear 
trend. Further simulations and future observations of clusters with temperature profiles and virial masses 



will be need to explore this issue further. 



4.2. The Isothermal /3-niodel 

The most often used x-ray mass estimator has been the isothermal /3-model which assumes that the 
gas is isothermal and that the gas density follows Equation |8[ Equation then becomes: 

M«.,.1.13xlO",j|^Jl(!i;l,Afe (.) 

assuming /z = 0.59. The /3-model mass can be written in terms of density contrast (Equation g) as: 

M{S, /3, TO = 1.1 X IQi'^r 1/2/33/2 (^^^ (^1 _ ^-01^) ''''^^ (l"^) 



We derive /3-model masses using the data of Fukazawa 1997 (hereafter F97) . In his study of the metal 



abundances and enrichment in the ICM, F97 presents a catalog of 38 clusters with temperatures, core 
radii, and f3 parameters derived from ASCA data (see Table g columns 2, 7, and 8). F97 estimated the 
P parameter from the ASCA GIS data using a Monte-Carlo method to take into account the spatial and 
energy dependence of the GIS PSF and estimated temperatures by excluding the central region of the x-ray 
emission to minimize cooling flow biases. 

The F97 data has the advantage of being a homogeneous sample, although ASCA is not the best 
instrument for surface brightness fitting due to its complicated PSF. However, the ASCA GIS /3-model fits 
are rather insensitive to the presence of cooling flows due to the high bandpass of the GIS. The clusters are 
also at low redshifts and so subtend a large area of the detector minimizing the effect of the poor resolution 
of the GIS. We compared the F97 (3 values with those of derived from ROSAT studies (primarily using the 



PSPC) of David et al. (1995), Cirimele et al. (1997)|, and various others taken from the literature via the 



compilation of [Arnaud fc Evrard (1998)| (see Figure |7| and Table g column 9). In general, they agree fairly 
well although ROSAT f3 values are higher by an average of about 5%. 

In Figure g, we plot the estimated /3-model mass at S — 200 (using Equation |o|) versus x-ray 
temperature. The relationship (M oc y^'i'Sio.os-j jg steeper than that seen using the EMN simulations, virial 
masses, or temperature profile masses. In addition, the relative normalization with respect to the other 
mass estimators is a function of density contrast. Increasing the density contrast shifts the /3-model masses 
lower with respect to the EMN (or virial mass) relation and closer to the temperature profile masses while 
decreasing the density contrast has the opposite effect. 

The isothermal /9-model implicitly assumes a dark matter density profile of p{r) ex r^^ in the outer parts 
of clusters. The density profiles of clusters in the EMN simulations and clusters with measured temperature 



profiles generally seem to follow a steeper profile p{r) oc r"^-^ (see Section 4.1 and Figure o). Therefore, the 
dark matter density profile implied by the isothermal /3 model is probably not an accurate representation of 
clusters, and, despite the agreement of the /3-profile with the observed gas density distribution, the different 
dark matter density profiles means that agreement between the isothermal /3 model estimate and other 
mass estimates will be a function of radius (hence density contrast). The /3-model masses underestimate 
the cluster mass at small radii and overestimate it at large radii. This has been pointed out previously by 
Bartelmann fc Steinmctz (1996)| and Markevitch et al. (1998)|. The different density profiles explain the 



relative steepness of the A/^ — Tx relationship. For an individual cluster, Afa/M(r(r)) oc r . At any given 



density contrast, r (x tJ and AI{T{r)) oc Tx , so the /3 model mass within that overdensity would be 
ex T^'^ similar to the fit given above. 



4-2.1. Gas Density Profiles 

Even if the underlying assumptions of the isothermal /3-model are flawed, Equation g| should still 
provide a good fit to the gas density distribution (see ^uto, Sasaki, fc Makino 1998 for a detailed discussion) . 



The EMN relation is basically a /3-model with f} = 0.68 for all clusters (depending on density contrast) 
regardless of x-ray temperature. However, the F97 data show a definite correlation of f3 with T^ (see 
Figure m. This indicates that the gas profile becomes shallower at lower masses. The relationship closely 



resembles that of the Metzler &: Evrard (1997) models with galactic winds, although with a normalization 
20% lower. 

The variation of (3 with T^ is unlikely to be an artifact of the F97 fitting procedure. As we noted 
in the previous section the F97 /3 values generally agree with ROSAT values. The correlation of /3 and 



Tx has also been noted previously in Einstein data David, Forman, fc Jones (1991) , and, more recently. 



Mohr fc Evrard (1997) have found a similar trend for (3 defined in a non-parametric and non-azimuthally 



averaged fashion using PSPC data. Arnaud & Evrard (1998) also note the behavior of /3 with T^ and the 



discrepancy between the /3-model and the expected EMN masses in their sample of clusters compiled from 
the literature. In fact, redoing the proceeding analysis with their sample gives virtually identical results. 



4.3. Surface Brightness Deprojection 

Another method for determining the mass of a cluster is x-ray image deprojection. The constraint of 
the observed surface brightness profile means that the profiles for the variables in Equation H & H can be 
determined by specifying one of them. The usual procedure is to divide the surface brightness emission into 
annuli. The outer pressure must be set in the outermost annuli (assumed to be due to gas not detected 
because its surface brightness is too low). The observed emissivity in the outer shell determines the 
temperature and hence the density. This procedure is then stepped inward and repeated. For more detailed 



discussion see Arnaud (1988), White et al. (hereafter WJF) and references therein. 



WJF present an analysis of 207 clusters using an x-ray image deprojection analysis of Einstein IPC 
and HRI data to estimate the masses of clusters. WJF choose the functional form of the gravitational 
potential as two isothermal spheres, representing the central galaxy and general cluster potentials. These are 
parameterized by a velocity dispersion and core radius. For the central galaxy, these are fixed at 350 km s~^ 
and 1 h^^ kpc. The velocity dispersion of the cluster potential is taken from the literature or interpolated 
from the x-ray temperature or luminosity using an empirical relation. The core radius is a free parameter in 
the analysis which, with the outer pressure, is constrained to produce a flat temperature profile. Therefore, 
the derived gravitational mass depends on the velocity dispersion, x-ray surface brightness distribution, and 
temperature. 

WJF determine the mass within the radius that they have x-ray data, which may be fairly small, 
while the EMN relation is only valid in the outer parts of the cluster (i.e. low values of density contrast). 
As with the G98 clusters, we can rescale the mass to a S = 200, at least for the WJF clusters which 
have data in the outer parts of the cluster. We have rescaled the clusters with data at density contrasts 



(= MwjF /{^/'■^T^'Ty^jpPc)) less than 2000 io 5 = 200 assuming p{r) oc r ^ as is the case for an isothermal 
sphere. 

In Figure M, we show the resulting M200 versus T^ relation. Although the lower Tx error bars are quite 
large, allowing the the points to be statistically consistent with the EMN relation, the relation is obviously 
steeper. Fitting this relation using the BCES method gives a M oc y^06±o.io somewhat steeper than than 
found from the F97 data. The scatter in mass around this relation is also larger, about 50%. 

Although WJF interpolated <j or T^ for many clusters, using only the 45 clusters with data at (5 < 2000 
which did not have interpolated data does not alter our results. Also, the results reported using a subsample 



of 19 clusters with better determined parameters by White & Fabian (1995) are consistent with the WJF 



results. Furthermore, a recent deprojection analysis by Peres et al. (1998) of 45 clusters with ROSAT PSPC 
or HRI data gives a similar relationship as the WJF data, although the Peres et al. data are generally at 
higher density contrasts (j> 3000) making the extrapolation to 5 — 200 even more uncertain. 

Interpretation of the deprojection results is difficult as the derived gravitational masses are a 
combination of optical (the velocity dispersion) and x-ray data (the core radius and temperature). However, 
like the /3-model masses the WJF masses follow a steeper relation due to the assumption of a potential 
which has a density profile which behaves as p{r) ex r^^. The large scatter is probably due to the use of the 
velocity dispersion to set the depth of the potential well. 



5. Conclusions 

We have examined the relationship between various galaxy cluster mass estimators and x-ray gas 
temperature. The resulting relationships generally agree to within <^ 40% in mass but with systematic 
offsets between different types of mass estimators. Using G98 virial masses and ASCA temperatures, we 
find good agreement with the EMN M-T relation after removing a few outliers. X-ray mass estimates using 
spatially resolved temperature profiles scale similarly but with a normalization about 40% lower. We note 



that the [Hjorth, Oukbir, fc Van Kampcn (1998)| lensing estimate mentioned at the end of Section g lies 



between these two estimates. 

Mass estimates based on the assumption of isothermality like the /3-model follow a steeper scaling 
relation due to the implicit assumption that the dark matter density scales p(r) oc r~^ at large radii while 
observational and numerical evidence suggests that clusters follow steeper profiles (i.e. p(r) oc r^^''). As a 
consequence, the /? model underestimates the mass at low radii and overestimates it at large radii. 

The intrinsic dispersion in the true M~T relationship is probably quite small. The scatter in the virial 
mass - temperature relation is consistent with most of the scatter being due to the dispersion in the virial 
mass estimator. The small scatter of the masses of clusters with spatially resolved temperature profiles also 
indicates that the dispersion in the M-T relation is probably <^ 10%. More insight may be gained through 
simulations of projection effects on various mass estimators, but the most comprehensive study of these 



effects Gen 1997 is limited to fairly low masses clusters (;^ 4 x lO^^/i ^Mq corresponding to Tx ^ 4 keV). 



Simulations covering a larger mass range are necessary to better constrain the scatter in the M-T relation. 

Our results also show that galaxy clusters seem to be affected by nongravitational processes such as 
energy injection by galactic winds. We find that (3, the asymptotic slope of the gas density profile, is a 
function of temperature. The dependence is consistent with models of energy injection by galactic winds 



which tend to make the gas profile shallower in less massive clusters Mohr fc Evrard (1997) . The steepening 



10 



of the cr-Tx relation is also consistent with energy injection by galactic winds (Bird et al. 1995). The good 

3/2 

agreement between the virial and temperature resolved mass estimates and a Tx scaling law indicates that 
any temperature changes (AT) due to energy injection must be small. However, meaningful constraints on 
AT are difficult given the relatively few clusters in this regime. 

In the future, more optical virial mass estimates for cooler clusters and larger samples of clusters with 
resolved temperature profiles will enable better constraints on the M-T relation. Better x-ray data will 
allow the effect of energy injection on the apparent temperature and the spatial distribution of the ICM 
to be disentangled. Future x-ray observations (e.g., with AXAF) will produce both the gas density and 
temperature maps of clusters and allow maps of the gas entropy to be constructed which should produce 
new constraints on energy injection histories and entropy variations within the cluster population. Such 
entropy maps can then be compared with those produced from simulations of various heating mechanisms. 

We would like to that D. White for provide machine readable tables from the WJF paper and Y. 
Fukazawa for providing a copy of his thesis. 
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Fig. 1. — Virial Mass- Temperature relation. The solid line is the theoretical EMN relation while the dashed 
line is a fit to the G98 virial masses (rescaled to an density contrast of 200) and ASCA temperatures. The 
error bars represent the 90% confidence intervals for temperature and 68% confidence for mass. The gray 
circles indicate clusters excluded from the fit (named in the plot). 
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Fig. 2. — Velocity Dispersion- Temperature relation. The solid line is cr oc Tx . The dashed line is a fit to 
the data excluding the clusters marked in gray. The error bars represent the 90% confidence intervals for 
temperature and 68% confidence for velocity dispersion. 
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Fig. 3. — Histogram of virial to best fit masses excluding outliers. 
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Fig. 4. — Ratio of virial to best fit masses as function of number of redshifts used to calculate virial mass. 
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Fig. 5. — Mass-Temperature Relation. Filled circles are clusters with masses measured using spatially 
resolved temperature profiles rescaled to an density contrast of 200. The dashed line is a fit to this data. 
Gray diamonds are isothermal /3-model masses within an overdensity of 200 for the F97 sample. The dotted 
line is the corresponding fit. The solid line is the theoretical EMN relation. Error bars are 90% confidence 
limits. 
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Fig. 6. — Density profiles of clusters with resolved temperature profiles as a function of radius, normalized 
to the EMN's r2oo- The solid line represents a fit to the data, S oc (r/r2oo)^^'**°'^- 
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Fig. 7. — Comparison F97 ASCA GIS P values to ROSAT /3 values. The circles are taken from David et al 
1995). The squares are from Cirimele et al. (1997), and the diamonds are taken from references listed in 



Arnaud fc Evrard (1998) 
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Fig. 8. — Comparison of viral and temperature resolved mass estimates. 
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Fig. 9. — Relation between (3 and T^ for F97 data. The dashed hne is a fit to the data (/3 ex T' 
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Fig. 10. — Mass-Temperature relation for WJF clusters using clusters with data for S < 2000. The solid 
circles are clusters which did not have any interpolated input data. The solid line is the theoretical EMN 
relation. The dashed line is a fit to the data. The error bars represent the 90% confidence intervals. 
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Table 1. Fitting Results M200 = cq x lO^^ y_ci /j-i^y^^ 



Sample 


Co 


Cl 


Comments 


EMN 


4.42 ±0.56 


1.5 




G98 


2.94 ±0.80 


1.71 ±0.20 


all clusters 


G98 


3.97 ±0.76 


1.53 ±0.13 


excluding outliers 


M{T^{r)) 


2.72 ±0.47 


1.48 ±0.12 




F97 


2.18 ±0.14 


1.78 ±0.05 




WJF 


2.02 ±0.29 


2.06 ±0.10 
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Table 2. Cluster Data 



Name 


Tx 


rVir 


Mvi. 


r-x 


Mx 


Tc 


/3 


Prosat 


Ref. 




keV 


h-^ Mpc 


IO^/i-^Mq 


h^^ Mpc 


IO^^/i-iMq 


ft-i Mpc 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


2A0335+096 


3.01±0.07 






0.50 


1.10 


0.023 


0.60 




3 


A0085 


6.31±0.25 


1.94 


9.88t?;g 






0.086 


0.60 


0.62±0.03 


10 


A0119 


5.59±0.27 


1.36 


2.50t°:?° 






0.231 


0.60 


0.56±0.04 


11 


A0194 


2.63±0.15 










0.069 


0.45 






A0262 


2.15±0.06 


1.05 


1 00+0.30 
-■-••^^-0.25 






0.032 


0.50 


0.53±0.03 


10 


A0399 


7.00±0.40 


2.23 


13.45t^;^^ 












1 


A0400 


2.31±0.14 


1.20 


2.49«:g 






0.051 


0.45 






A0401 


8.00±0.40 


2.30 


13.69+^-^« 












1 


A0426 


6.79±0.12 


2.05 


9.08t?:Ji 






0.020 


0.45 






A0478 


6.90±0.35 










0.077 


0.70 






A0496 


4.13±0.08 


1.37 




0.50 


1.50 


0.035 


0.55 




4 


A0539 


3.24±0.09 


1.26 


r, Q1+0.52 






0.082 


0.60 


0.65±0.03 


10 


A0754 


9.50±0.60 


1.32 


4 OQ + 1.04 
^•^•J-0.72 












1 


A1060 


3.24±0.06 


1.22 


1 qn+0-38 


1.00 


1.80 


0.040 


0.55 


0.61 


5,9 


A1656 


8.38±0.34 




. qy+0.68 






0.208 


0.75 


0.75±0.03 


9 


A 1795 


5.88±0.14 


1.67 


5.86+;-^^ 


1.00 


3.00 


0.068 


0.65 


0.74±0.02 


4,10 


A2029 


9.10±1.00 


2.33 


6.82t-.« 


0.96 


4.71 








1,7 


A2063 


3.68±0.11 


1.33 


3.04«:f3 






0.067 


0.60 


0.67±0.04 


10 


A2107 


3.78±0.19 


1.24 


2.62«:«4 












2 


A2142 


9.70±1.30 


2.26 


17 8/1+3-99 
l'-84_3.49 












1 


A2147 


4.91±0.28 










0.109 


0.50 






A2163 


13.83±0.76 






1.00 


10.70 
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A2199 


4.10±0.08 


1.60 


5.7ltl-J? 


0.50 


1.28 


0.040 


0.60 


0.62±0.02 


4,11 


A2256 


7.08±0.23 


2.70 


23.12+^-^^ 


1.50 


6.00 


0.228 


0.75 


0.81±0.01 


6,9 


A2319 


8.90±0.34 


3.09 


39.541^:85 






0.135 


0.55 






A2634 


3.70±0.28 


1.40 


A Q1+1-35 






0.123 


0.50 


0.58±0.02 


11 


A2670 


4.45±0.20 


1.70 


c cc+0.72 
0-00-0.58 












2 


A3266 


8.00±0.50 


2.21 


11.7011:^^ 












1 


A3376 


4.00±0.40 


1.38 


3.64t°-i« 












1 


A3391 


5.40±0.60 


1.33 


o f;i+2.20 
•J-OJ^-1.35 












1 


A3395 


5.00±0.30 


1.70 


b.mtXil 












1 


A3526 


3.68±0.06 






0.50 


0.80 


0.038 


0.50 




3 


A3558 


5.12±0.20 


1.95 


ii-54li:»? 






0.075 


0.50 


0.61 


9 


A3571 


6.73±0.17 


2.09 


8.17t^;f9 






0.086 


0.60 






A3667 


7.00±0.60 


1.94 


ii.75tl:«° 












1 


A4059 


3.97±0.12 






0.50 


1.50 


0.075 


0.65 




3 


AWM4 


2.38±0.17 










0.034 


0.55 


0.47 




AWM7 


3.75±0.09 


1.73 


5.77tl-J« 


1.00 


2.00 


0.062 


0.55 


0.53±0.01 


5,9 


HCG51 


1.23±0.06 










0.058 


0.55 






HCG62 


1.05±0.02 










0.003 


0.40 






HYDRA-A 


3.57±0.10 










0.036 


0.60 






MKW3S 


3.68±0.09 






1.00 


2.00 


0.047 


0.65 




3 


MKW4 


1.71±0.09 


1.05 


1 -1^+0.35 
-'-•-'-O-0.26 






0.009 


0.45 






MKW4S 


1.95±0.17 










0.023 


0.40 






MKW9 


2.23±0.13 










0.025 


0.45 
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Table 2 — Continued 



Name 


Tx 


■^vir 


Mvi. 


^x 


Mx 


fc 


/3 


Prosat 


Ref. 




keV 


ft-i Mpc 


10"/i-^Mq 


ft-i Mpc 


IO^^/i^^Mq h-^ Mpc 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NGC2300 


0.88±0.03 












0.024 


0.40 






NGC5044 


1.07±0.01 












0.009 


0.50 






NGC507 


1.26±0.07 












0.014 


0.45 






OPHIUCHUS 


10.26±0.32 












0.113 


0.60 






S753 


2.23±0.12 


1.07 


1 01+0.67 
J^-jJ^-0.50 














2 


TRIAUST 


10.05±0.69 












0.126 


0.60 


0.63±0.02 


9 


VIRGO 


2.58±0.03 


1.26 


9 (]A+0-'^» 








0.007 


0.40 


0.46 


9 



References. — (1) Tx from Markevitch (1998) (2) Tx derived for tliis paper (3) Mx & r^ from Ohashi (1997) (4) Mx & rx 
from Mushotzky ot al. (1994) (5) M^ & r^ from Loewenstein & Musliotzky (1996) (6) M^ & r^ from Markevitch & Vikhlinin 
(1997) (7) Mx & rx from Sarazin, Wise, & Markevitch (1998), their value of T^ = 8.691!^;^^ is used in Section 4.1 (8) M^ & r^ 
from Markevitch et al. (1996) (9) Prosat from Arnaud & Evrard (1998) (10) Prosat from David, Jones, & Forman (1995) 
(11) Prosat from Cirimele, Nesci, & Trevesc (1997) 



